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ABSTRACT 

We analyse two ASCA observations of the highly eccentric 09 III + Bl III binary 
Iota Orionis obtained at periastron and apastron. Based on the assumption of a strong 
colliding winds shock between the stellar components, we expected to see significant 
variation in the X-ray emission between these phases. The observations proved other- 
wise: the X-ray luminosities and spectral distributions were remarkably similar. The 
only noteworthy feature in the X-ray data was the hint of a proximity effect during 
periastron passage. Although this 'flare' is of relatively low significance it is supported 
by the notable proximity effects seen in the optical (Marchenko et al. 2000) and the 
phasing of the X-ray and optical events is in very good agreement. However, other 
interpretations are also possible. 

In view of the degradation of the SIS instrument and source contamination in the 
GIS data we discuss the accuracy of these results, and also analyse archival ROSAT 
observations. We investigate why we do not see a clear colliding winds signature. A 
simple model shows that the wind attenuation to the expected position of the shock 
apex is negligible throughout the orbit, which poses the puzzling question of why the 
expected 1/ D variation (i.e. a factor of 7.5) in the intrinsic luminosity is not seen in the 
data. Two scenarios are proposed: either the colliding winds emission is unexpectedly 
weak such that intrinsic shocks in the winds dominate the emission, or, alternatively, 
that the emission observed is colliding winds emission but in a more complex form 
than we would naively expect. Complex hydrodynamical models are then analyzed. 
Despite strongly phase-variable emission from the models, both were consistent with 
the observations. We find that if the mass-loss rates of the stars are low then intrinsic 
wind shocks could dominate the emission. However, when we assume higher mass- 
loss rates of the stars, we find that the observed emission could also be consistent 
with a purely colliding winds origin. A summary of the strengths and weaknesses of 
each interpretation is presented. To distinguish between the different models X-ray 
observations with improved phase coverage will be necessary. 



Key words: stars: individual: Iota Orionis (HD 37043) 
general - X-rays: stars 



■ stars: early-type - binaries: 



1 INTRODUCTION 



* Charge de Recherches FNRS, Belgium 
© 2000 RAS 



The Orion OBI stellar association is one of the brightest 
and richest concentrations of early-type stars in the vicinity 
of the Sun (Warren & Hesser 1977) . It contains large num- 
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bers of young O- and B-type stars, and this, combined with 
its position well below the galactic plane (< 6 > = -16°) 
with subsequent low foreground absorption, makes it an op- 
timal object to study. The OBI association is also famous 
for an exceptionally dense concentration of stars known as 
the Trapezium cluster, located near 9^ Ori. The stellar den- 
sity of the Trapezium cluster has been estimated by several 
authors in recent years and its central region is now thought 
to exceed ~ 4 x 10'' stars pc~^ (McCaughrean & Stauffer 
1994), making it one of the densest young clusters currently 
known. 

Iota Orionis (HR 1889; HD 37043) is a well known 
highly eccentric (e = 0.764) early-type binary system (09 III 
+ Bl III, P = 29.13376 d) located in the C4 subgroup of the 
Orion OBI association, approximately 30 arcminutes South 
of the Trapezium cluster. In the last decade Iota Orionis has 
drawn attention for the possibility of an enhanced, focused 
wind between the two stars during the relatively close perias- 
tron encounter (Stevens 1988; Gies, Wiggs & Bagnuolo 1993; 
Gies et al. 1996). Although data consistent with this effect 
was presented (residual blue-shifted emission in Ha that ac- 
celerates from about -100 kms"'^ to -180 kms"'^ ), previous 
evidence was inconclusive, and the possibility of a random 
fluctuation of the primary's wind was not ruled out by Gies 
et al. (1996). These authors also found no evidence of any 
systematic profile deviations that are associated with non- 
radial pulsations. The latest optical monitoring (Marchenko 
et al. 2000), again obtained data consistent with a proximity 
effect. 

Iota Orionis is perhaps even more interesting for the 
possible existence of a strong colliding winds interaction re- 
gion between the two stars, which should be most clearly 
recognized at X-ray wavelengths. As demonstrated by Cor- 
coran (1996), perhaps the best direct test for the impor- 
tance of colliding winds emission is to look for the expected 
variation with orbital phase. X-ray data can also be used 
to extract information on the characteristics of this region 
(e.g. geometry, size, temperature distribution), which can 
in turn constrain various stellar parameters including mass- 
loss rates. This possibility has already been explored for the 
Wolf-Rayet (WR) binary 7^ Velorum (Stevens et al. 1996). 
Furthermore, a direct comparison of the recently presented 
sudden radiative braking theory of Gayley et al. (1997) with 
observed X-ray fluxes has not yet been made. For this goal. 
Iota Orionis has the dual advantage of a high X-ray flux 
(partially due to its relative proximity, D ~ 450 pc) and 
a highly eccentric orbit (the varying distance between the 
stars should act as a 'probe' for the strength of the radia- 
tive braking effect) . Iota Orionis is additionally intriguing for 
its apparent lack of X-ray variability (see Section H) whilst 
colliding wind theory predicts a strongly varying X-ray flux 
with orbital phase. 

In this paper we present the analysis and interpreta- 
tion of two ASCA observations of Iota Orionis proposed by 
the XMEGA group, as well as a reanalysis of some archival 
ROSAT data. These are the first detailed ASCA observa- 
tions centered on this object, although ASCA has previ- 
ously observed the nearby O-stars S Ori and A Ori (Corco- 
ran et al. 1994) and surveyed the Orion Nebula (Yamauchi 
et al. 1996). ASCA observations of WR-l-0 binaries (e.g. 
7^ Velorum - Stevens et al. 1996; WR 140 - Koyama et 
al. 1994) have provided crucial information of colliding stel- 



lar winds in these systems. The new ASCA observations 
of Iota Orionis enable us to perform a detailed study of 
the X-ray emission in this system. The periastron obser- 
vation covered ~ 65 ksec and includes primary minimum 
{(j) = —0.009), periastron passage {(j> — 0.0), and quadra- 
ture (0 — 0.014), and allows an unprecedented look at in- 
teracting winds in an eccentric binary system over a range 
of different orientations. The apastron observation covered 
~ 25 ksec and allows us to compare the X-ray properties at 
the maximum orbital separation. Thus these two new ob- 
servations offer the additional benefit of improving the poor 
phase sampling of this object. 

The new ASCA observations were coordinated with an 
extensive set of ground based optical observations reported 
in Marchenko et al. (2000). In this work, the spectra of the 
components were successfully separated, allowing the refine- 
ment of the orbital elements (including the restriction of the 
orbital inclination to 50° ^ i ;$ 70°) and confirmation of 
the rapid apsidal motion. Strong tidal interactions between 
the components during periastron passage and phase-locked 
variability of the secondary's spectrum were also seen. How- 
ever, no unambiguous signs of the bow shock crashing onto 
the surface of the secondary were found. These results are 
extremely relevant to the interpretation of the ASCA X-ray 
data. 

This paper is organized as follows: in Section m we 
briefiy discuss previous X-ray observations of Iota Orionis; in 
Sections we present our analysis of the new ASCA datasets 
and reexamine archival .ROSAT data. In Section M we exam- 
ine possible interpretations of the data. Comparisons with 
other colliding wind systems are made in Section 0, and in 
Section we summarize and conclude. 



2 PREVIOUS X-RAY OBSERVATIONS OF 
IOTA ORIONIS 

On account of its high X-ray flux. Iota Orionis has been ob- 
served on several occasions. X-ray emission from Iota Orio- 
nis was reported by Long & White (1980) in an analysis of 
Emstem Imaging Proportional Counter (IPC) data, and a 
count rate of 0.289 ± 0.007 cts/s in the 0.15-4.5 keV band 
was obtained. The spectral resolution (E/AE) of the IPC 
was only 1-2, so only the crudest spectral information could 
be obtained, although it was observed that the spectrum was 
soft, peaking below 1 keV. The authors noted that qualita- 
tively good fits could be obtained with a variety of spectral 
models leading to an uncertainty of a factor of 2 in their 
quoted luminosity of 2.3 X 10^^ ergs"^ (0.15-4.5 keV). 

Since this observation the X-ray properties of Iota Ori- 
onis have been studied on numerous occasions (Snow, Cash 
& Grady 1981; CoUura et al. 1989; Chlebowski, Harnden 
& Sciortino 1989; Waldron 1991; Gagne & Caillauff 1994; 
Geier, Wendker & Wisotzki 1995; Berghofer & Schmitt 
1995a; Kudritzki et al. 1996; Feldmeier et al. 1997a). Due 
to the different X-ray satellites, data-analysis techniques, 
and orbital phases, meaningful comparisons between the 
datasets are difficult, and we simply refer the reader to the 
above mentioned papers. We note, however, the results of 
previous X-ray variability studies of Iota Orionis, performed 
by Snow et al. (1981) and CoUura et al. (1989) on Einstein 
data, and by Berghofer & Schmitt (1995a) on ROSAT datin. 
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Table 1. Einstein IPC observations of Iota Orionis. 



Sequence 


Date of 


Phase 




IPC counts/s 




number 


obs. 





(LW80) 


(SCG81) 


(C89) 


3842 


23/9/79 


0.34 


0.289 






11263 


03/3/80 


0.90 




0.609 




5095 


05/3/80 


0.97 




0.585 


0.417 


5096 


23/3/80 


0.59 




0.607 


0.413 


10413 


16/2/81 


0.92 






0.367 



Notes. The count rates hsted were obtained from the follow- 
ing papers: LW80 - Long & White (1980); SCG81 - Snow, Cash 
feCrady (1981); C89 - CoUura et al. (1989). 

Snow et al. found that the count rates from 3 IPC observa- 
tions were twice as high as the Long & White observation, 
leading to the simple interpretation that the X-ray emission 
from Iota Orionis could vary by up to a factor of two (with 
a statistical significance of 3cr), though it was noted that at 
the time of the observations there were problems with the 
instrumental gain shifts of the IPC. CoUura et al. (1989) 
reexamined the issue of variability with a more rigorous 
method and found that they could only place upper lim- 
its of 12% on the amplitude of variability over timescales 
ranging from 200 s up to the duration of each observational 
interval. In contrast to the findings of Snow et al. (1981) no 
evidence of long-term variability between observations was 
detected. Finally Berghofer & Schmitt (1995a) found an es- 
sentially flat lightcurve from ROSAT observations, consis- 
tent with the previous findings of Collura et al. (1989), al- 
though again the phase-sampling was rather poor. The lack 
of dramatic variability in these previous studies could be a 
result of poor phase sampling, and/or the soft response of 
the instruments (assuming that X-rays from the wind col- 
lision are harder than those intrinsically emitted from the 
individual winds). Our ASCA observations were specifically 
designed to address this by i) sampling at phases where the 
variation should be maximized; and ii) by leveraging the 
extended bandpass of ASCA to measure the hard flux. 

Table hi summarizes the Einstein IPC observations and 
the rather confusing situation with regard to whether and at 
what level there is any variability (note that it is not meant 
to imply that there is any real source variability). It also 
emphasizes the poor phase sampling and serves as a useful 
highlight as to one of the reasons why the latest ASCA data 
presented in this paper were obtained. In column 3 we list 
the orbital phase at the time of observation, calculated from 
the ephemeris 

HJD {periastron) = 2451121.658 + 29.13376B 

derived from the latest optical monitoring of this system 
(Marchenko et al. 2000). As a result of the large number of 
orbits that have elapsed the listed phases are only approxi- 
mate and do not account for the apsidal motion detected in 
Iota Orionis (Stickland et al. 1987, Marchenko et al. 2000) . 
We also draw attention to the different count rates obtained 
by different authors when analysing the same datasets, and 
the relatively constant count rates obtained from different 
datasets by the same authors. 

A further ROSAT analysis of Iota Orionis was published 
by Geier et al. (1995), who obtained a position-sensitive 



proportional counter (PSPC) observation centered on the 
Trapezium Cluster. This discovered that most of the X-ray 
emission from the region originated from discrete sources, 
in contrast to the previous Einstein data where the sources 
were spatially unresolved. Nearly all of the 171 X-ray sources 
were identifled with pre-main sequence stars in subgroups 
Ic and Id of the Orion OBI association. The 4 point- 
ings made from March 14-18 1991, together gave a total 
effective exposure time of 9685 s. However, Iota Orionis 
was significantly off-axis in all of these pointings. Assum- 
ing a Raymond & Smith (1977) thermal plasma spectral 
model, a single-temperature fit gave kT = 0.22 keV and 
Nh ~ lO^^cm"^, whilst a two-temperature fit with a single 
absorption column gave kTi = 0.12 keV, kT2 — 0.82 keV, 
and Nh ^3x 10^°cm~^ 

Two analyses of ROSAT PSPC spectra which ignored 
the binary nature of Iota Orionis have recently been pub- 
lished by Kudritzki et al. (1996) and Feldmeier et al. (1997a). 
The X-rays were assumed to originate from cooling zones be- 
hind shock fronts. These are a natural development of the 
intrinsically unstable radiative driving of hot-star winds (see 
e.g. Owocki, Castor & Rybicki 1988). Although the authors 
stated that good fits were obtained, no formal indication of 
their goodness was given. Because binarity may significantly 
affect the X-ray emission, we use this as the basis of our in- 
terpretation of the data. This approach complements these 
alternative single-star interpretations. 



3 THE X-RAY DATA ANALYSIS 

Despite previous attempts to characterize the X-ray proper- 
ties of Iota Orionis, our understanding remains poor. For in- 
stance, from the previous Einstein and ROSAT data, it looks 
like a typical O-star. However, on account of its stellar and 
binary properties. Iota Orionis should reveal a clear X-ray 
signature of colliding stellar winds which one would expect 
to show dramatic orbital variability as the separation be- 
tween the stars changes. ASCA should be well suited for the 
study of this emission on account of its high spectral resolu- 
tion and greater bandpass. In this section we report on the 
analysis performed on the two new phase-constrained ASCA 
datasets. This is followed by a reanalysis of two archival 
ROSAT datasets, and a comparison with results obtained 
from colliding wind models. 



3.1 The ASCA analysis 

3.1.1 Data reduction 

Iota Orionis was twice observed with ASCA during 19971']. 
An observation on 21 September 1997 was timed to coin- 
cide with periastron passage, whilst the 6 October observa- 
tion was timed to coincide with apastron passage. The event 
files from each instrument were screened using the FTOOL 
ASCASCREEN. For the standard SIS analysis the BRICHT 
datamode was processed using medium and high bit rate 
data. Hot and fiickering pixels were removed and the stan- 
dard screening criteria applied. 

T Details of the ASCA satellite may be found in Tanaka, Inoue 
& Holt (1994). 
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Iota Orionis Periastron Observation — E < 1.0 keV 



Iota Orionis Periastron Observation — E > 1.0 keV 
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Figure 1. Screened images from the SISO periastron data at energies below 1.0 keV (left) and above 1.0 keV (right). The serendipitous 
source is noticeably brighter in the image on the right, which tells us that its spectrum is relatively harder than Iota Orionis. The left 
image clearly shows the characteristic 'Maltese Cross' of the XRT PSF. 



Table 2. Effective exposure times and number of counts (back- 
ground corrected) for each of the four instruments onboard ASCA 
during the periastron and apastron pointings. The percentage of 
counts in the source region due to background varied between 20- 
25 per cent for the periastron pointing and 15-23 per cent for the 
apastron pointing. The orbital phase spanning each observation 
is also listed in column 1. 



Observation 


Inst. 


Exp. 


Cts in 


Ct rate (s'^) 


(Phase, <j)) 




(s) 


source 
region 


in source 
region 


Periastron 


SISO 


65672 


12543 


0.191 


0.9905-1.068 


SISl 


66545 


9729 


0.146 




GIS2 


72288 


4301 


0.060 




GIS3 


72336 


5085 


0.070 


Apastron 


SISO 


23934 


4856 


0.203 


1.502 - 1.529 


SISl 


23798 


3225 


0.136 




GIS2 


24496 


1416 


0.058 




GIS3 


24496 


1814 


0.074 



The four X-ray telescopes (XRTs) onboard ASCA each 
have spatial resolutions of 2.9 arcmin half power diameter. 
In spite of the broad point-spread function (PSF), the jit- 
tering of the spacecraft can appear on arcminute scales. To 
avoid inaccurate flux determinations and spurious variabil- 
ity in the hght curves, the radius of the source region should 
normally be no smaller than 3 arcmin [ASCA Data Reduc- 
tion Guide). For most cases the recommended region filter 
radius for bright sources in the SIS is 4 arcmin. Because of 
concerns of contamination from nearby sources, our stan- 
dard analysis adopted a source radius of 3 arcmin. 

The estimation of the X-ray background initially proved 



somewhat troublesome. The presence of a nearby source 
forced a thin annulus, creating a large variance on the back- 
ground spectrum, which invariably still contained contam- 
inating counts. The available blank sky backgrounds were 
also found to be unsatisfactory. After much consideration, 
the background was estimated from the entire CCD exclud- 
ing the source region and other areas of high count rate. This 
method produced a large number of background counts, re- 
ducing the uncertainty in its spectral shape, and resulting 
in a spectral distribution for the source which was better 
constrained. 

For the GIS analysis, the standard rejection criteria for 
the particle background were used together with the rise 
time rejection procedure. Unlike the SIS, the intrinsic PSF's 
of the GIS are not negligible compared to the XRT's and 
a source region filter of 4 arcmin radius was used. The 
serendipitous source detected by the SIS instruments was 
not spatially resolved by the GIS instruments and hence the 
source spectra extracted from the latter contain contami- 
nating counts from this object. 

Two methods of background subtraction for the GIS in- 
struments are commonly used. Either one can use blanksky 
images with the same region filter as used for the source ex- 
traction, or one can choose a source-free area on the detector 
at approximately the same off-axis angle as the source. In 
our analysis the second method was favoured due to the 
following reasons: 

(i) The cut-off-rigidity (COR) time-dependence can be 
correctly taken into account. 

(ii) The blanksky background files were taken during the 
early stages of the ASCA mission, and do not include the 
secular increase of the GIS internal background. 

(iii) The blanksky background files were also taken from 
high Galactic latitude observations, and hence possible dif- 
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ROSAT PSPC Observation of Iota Orionis 
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Figure 2. flOSylT PSPC image of the Iota Orionis field. Tliis 
image number is rf200700n00, and clearly shows the X-ray sources 
near Iota Orionis, which is the brightest source at the centre of the 
field. The serendipitous ASCA source is labelled. Contour levels 
are 1, 2, 5, 10, 20 and 50 counts. Epoch 2000.0 coordinates are 
shown. 



the background-subtracted source was constant gave Xv = 
0.85. For the apastron observation a corresponding analysis 
gave xt = 0.75. If the majority of the X-ray emission was 
from the wind coUision region, one would expect significant 
variation given the highly eccentric nature of the system. 
There is also little difference between the count rates of the 
periastron and apastron observations. 

However, close to periastron passage (corresponding to 
2.38 X 10* s after the start of the observation) there does 
seem to be a small transient spike in both the SISO (see 
Fig. |3|) and GIS2/3 count rates (although nothing appears 
in the SISl lightcurve). This 'flare' is of relatively low signif- 
icance, and could be entirely instrument related as it occurs 
just before an interval of bad data. Supporting evidence for 
this comes from Moreno & Koenigsberger (1999) who have 
recently investigated the effect of tidal interactions between 
the stars around periastron passage. They found that if the 
radius of the primary is smaller than 15 R©, the enhance- 
ment in the mass-loss rate at periastron, and therefore the 
effect on the colliding winds emission, should not be signifi- 
cant. However, the phasing of the X-ray fiare is almost per- 
fect with a similar event seen in the optical (Marchenko et 
al. 2000). To examine whether the fiare possibly occurred in 
the background a cross-correlation analysis was performed. 
No significant correlation was found between the source and 
background lightcurves, whilst a strong correlation clearly 
occurred between the unsubtracted and subtracted source 
lightcurves. We conclude, therefore, that it is unlikely that 
this 'fiare' occurs in the background data. 



fuse X-ray emission near Iota Orionis cannot be taken into 
account. 

Despite these reasons, the blanksky method of background 
subtraction was also tested and in contrast to the SIS analy- 
sis both methods produced spectral fit results which agreed 
within their uncertainties. However, all subsequent work in- 
cluding the results reported in this paper used a background 
extracted from the same field as the source. 

The effective exposure times and count rates for the var- 
ious instruments during the periastron and apastron point- 
ings are detailed in Table Q. The numbers in columns 4 and 
5 refer only to the specific source regions extracted. 

In Fig. nl we show the SISO field from the periastron 
observation for two different energy bands. In addition to 
Iota Orionis, a serendipitous source is clearly visible to the 
right, being most noticeable in the harder image. A ROSAT 
PSPC image of Iota Orionis is shown in Fig. H, and demon- 
strates more clearly the crowding of sources in this region 
of the sky. The serendipitous source in Fig. hi is tentatively 
identified as NSV 2321, an early G-type star (optical coords: 
RA(2000.0) = 05'' 35™ 22=, Dec.(2000.0) = -05° 54' 35")- 
As already noted, this source complic ates t he analysis, of 
which more details are given in Section 3.1.4. 



3.1.2 The X-ray lightcurves 

Lightcurves from all 4 instruments were extracted for each 
of the two observations and analysed with the XRONOS 
package. A subset of these is shown in Fig. H. Somewhat 
surprisingly the lightcurves were found to be remarkably 
constant. For the periastron SISO observation, a fit assuming 



3.1.3 The X-ray spectra 

Using the XSPEC package, source X-ray spectra were ex- 
tracted from each dataset, re-binned to have a minimum 
of 10 counts per bin (as required for x^-fitting), and fitted 
with both single and two-temperature Raymond-Smith (RS) 
spectral models. We emphasize at this point that this sim- 
ple analysis is an effort to characterize i) the overall shape 
of the spectrum and its temperature distribution, ii) the 
amount of variability in the spectral parameters, and iii) 
the discrepancies between theoretical colliding stellar wind 
spectra and the observed data (see Section 4.2). It is not 
meant to imply that in the two-temperature fits the emis- 
sion physically occurs from two distinct regions at separate 
temperatures, and is simply in keeping with X-ray analy- 
sis techniques commonly used today. The upper panel of 
Fig. H shows a two-temperature RS model with one absorp- 
tion component that has been simultaneously fitted to the 
SISO and SISl periastron datasets, whilst the lower panel 
shows the corresponding fit to the GIS2 and GIS3 perias- 
tron datasets. We include data down to ~ 0.5 keV for the 
SIS instruments (important for measuring changes in Nh, 
although see Section 3.1.4) and to ~ 0.8 keV for the GIS 
instruments (where the effective area is roughly 10 per cent 
of the maximum effective area for this instrument, thus en- 
abling us to extract maximum information about the emis- 
sion below 1.0 keV). 

Some line emission is clearly seen, in particular at 
~ 1.85 keV and at ~ 2.30 keV, of which the former is most 
probably from SiXIII. Numerous spectral models were fit- 
ted to the data on both an instrument by instrument basis 
and to various combinations of the four instruments includ- 
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Figure 3. Background subtracted lightcurves from the SISO detector for the periastron (top) and apastron (bottom) observations. They 
are remarkably constant considering the anticipated wind interaction. However, an emission 'spike' is clearly seen in the periastron 
lightcurve, which occurs just before the point of closest approach of the stars (MJD = 50713.28536). It is not seen in the SISl data, 
though, and we note that it occurs just before a period of bad data, and therefore could be entirely instrumental. 
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ing all four instruments together. First we attempted sin- 
gle temperature solar abundance fits. Only the fit to the 
combined GIS2 and GIS3 apastron data was acceptable: 
Nh = 2.71 ± 0.63 X 10^^ cm-2, kT = 0.61 ± 0.06 keV, 
EM = 7.7 ± 1.7 X 10^* cm"^, xl = 0.92. Fits using non- 
solar abundances were not significantly improved, although 
the addition of a second independent plasma component did 
result in much better fits and we report some results in Ta- 
ble H. Immediately obvious is that the fits to the periastron 
and apastron SIS data are very similar, both in terms of the 
spectral parameters and the resulting luminosities. Clearly 
there is not the order of magnitude variation which we ex- 
pected in the intrinsic luminosity. This is also true for the 
corresponding fits to the GIS data. Thus it appears at first 
glance that any variability is small. This has implications 
for a colliding winds interpretation of the data, although 
at this stage we do not rule this model out. Other general 
comments which we can make on the spectral fits are: 

(i) The SIS datasets return much higher values of Nh 
and much lower values of the characteristic temperature, 
kT, than the GIS datasets. 

(ii) Both the SIS and GIS fits return absorbing columns 
which are greater than the estimated interstellar column 
(2.0 X 10^° cm"^; ShuU & Van Steenberg 1985, and Savage 
et al. 1977) arguing for the presence of substantial circum- 
stellar absorption consistent with strong stellar winds. 

(iii) Due to the larger variance on the GIS data, the GIS 
fit results are statistically acceptable, whereas those to the 
SIS data are not. 

(iv) If the global abundance is allowed to vary during the 
fitting process, the general trend in the subsequent fit results 
is that the characteristic temperature increases, the absorb- 
ing column decreases, the normalization increases, and the 
abundance fits at 0.05-0.2 solar. The resulting xi is basi- 
cally unchanged. These fit results bear all the hallmarks of 
the problems mentioned by Strickland & Stevens (1998) and 
we question their accuracy. 

(v) The metallicities generally fit closer to solar values if 
a two-temperature RS model is used. 

(vi) If a simultaneous SIS and GIS fit is made, the re- 
turned temperatures are similar to those from the SIS 
dataset alone. This is perhaps not so surprising, however, 
given that the variance on the GIS data is much larger than 
the SIS, and therefore that the SIS data has a much larger 
infiuence in constraining the fit. 

Two-temperature RS models with separate absorbing 
columns to each component were also fitted, as were two 
temperature spectral models with neutral absorption fixed 
at the ISM value and an additional photoionized component. 
However, both failed to significantly improve the fitting and 
we do not comment on these further. 



3.1.4 Problems with the analysis 

Although the differences between the individual SISO and 
SISl datasets, and between the GIS2 and GIS3 datasets are 
within the fit uncertainties, there is a large discrepancy be- 
tween the fit results made to the SIS datasets on the one 
hand and the GIS datasets on the other. Models which fit 
the GIS do not fit the SIS, and vice- versa. Whilst we can 
already state with some confidence that the X-ray emission 
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Figure 4. Two-temperature Raymond-Smith spectral model fits 
to the combined SISO and SISl data (upper panel) and the com- 
bined GIS2 and GIS3 data (lower panel) from the entire peri- 
astron observation. Both spectral components and the combined 
model are shown, together with the data. A minimum of 10 counts 
is in each data bin. 



is not strongly variable, the above is obviously a large cause 
for concern which we should still explore. For instance, if 
we ignore the lack of variability for a moment and assume 
that a colliding winds scenario is the correct model, then 
the fitted temperatures should provide some information on 
the pre-shock velocities of the two winds, and whether these 
vary as a function of orbital phase. Hence for this reason we 
wish to determine which of the SIS or GIS temperatures is 
most accurate. 

First the fit-statistic parameter space was investigated 
for false or unphysical minima. This is shown in Figure n 
for a single-temperature RS model fitted to the SISO and 
GIS2 periastron data. It is clear that there are no additional 
minima between the two shown, and the SISO minimum in 
particular is very compact. However, at very high confidence 
levels there is a noticeable extension of the GIS2 confidence 
region towards the SISO minimum. Despite this the results 
seem to be totally incompatible with each other. Other pos- 
sible explanations for the observed discrepancy are offered 
below: 

Degradation of the SIS instruments 
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Table 3. Two-temperature spectral fitting results. Fits made to the combined SISO and SISl datasets and the combined GIS2 and GIS3 
datasets are shown, as well as that to an archive ROSAT dataset (see Table U). The estimated ISM column is 2.0 X 10 cm~ . The 
emission measure and luminosities (0.5 — 2.5 keV) were calculated assuming a distance of 450 pc. The emission measure to the soft 
component of the rf200700n00 fit could not be constrained and led to unrealistic luminosities. Consequently we do not include it here. 
Many other spectral models were also examined (see main text) but for conciseness and simplicity wc do not report them here. 



Data 



Nh 
(lO^icm 



fcTi 
(keV) 



kT2 
(keV) 



(lO^^cm- 



EM2 
(10^4 cm-3) 



(DOF) 



La;(int) 
(lO^^ergs-i ) 



La;(abs) 
(10^2 erg s-i) 



Periastron SIS 
Apastron SIS 
Periastron GIS 
Apastron GIS 
rp200700n00 



5 07+0-28 
3-^'_0.28 
•x qn+0'53 
o.»u_g j3 

7+0.55 

-0.76 

4 73+1.43 

^•''^-1.46 

j+0.42 

-0.15 



1.97^ 



0.28^ 



o.i4t^;^l 


nfii+002 


260+^} 


40+0.5 


2.06 (225) 


20.17 


1.04 


o.i5t^;Sl 


nfii+003 


94+47 
y*-47 


4 1+0-7 
^■-'--0.7 


2.06 (143) 


8.99 


1.00 


61+°"^ 

'-'■"1-0.04 


o 1 9+2.31 
■J-l^_0.92 


6.06+?-^I 


94+0-36 


1.08 (381) 


2.01 


1.08 


9c;+oio 


2.i5tS;:I^ 


68t^f 


9 n+0.96 


0.77 (210) 


12.3 


1.46 


11+003 


o.7ot«:li 


26.7t?« 


ogn+0.73 
■'•"'^-0.73 


0.61 (11) 


1.85 


1.60 



The HEASARC website (http://heasarc.gsfc.nasa.gov/ 
docs/asca/ watchout.html) reports that there is clear evi- 
dence for a substantial divergence of the SISO and SISl 
detectors since late 1994, with even earlier divergence be- 
tween the SIS and GIS data. It is also noted that both the 
SISO and SISl efficiencies below 1 keV have been steadily 
decreasing over time, which at 0.6 keV can be as much as 20 
per cent for data taken in 1998. The loss in low-energy ef- 
ficiency manifests itself as a higher inferred column density, 
introducing an additional Nh uncertainty of a few times 
10^" cm~^. This is consistent with our results where the SIS 
detectors return higher values of Nh than the GIS (see Ta- 
ble H), although we note that the difference in our fits can 
often be much larger. An increasing divergence of the SIS 
and GIS spectra in the energy range below 1 keV has also 
been noted by Hwang et al. (1999). As also noted on the 
HEASARC website a related problem is that the released 
blank-sky data cannot be used for background-subtraction 
for later SIS data since the blank-sky data was taken early 
on in the ASCA mission. This is als o consistent with the 
problems mentioned in Section p.l.l| in obtaining good fits 
when using this method of background subtraction. Whilst 
these problems offer a convincing explanation for the dis- 
crepancies in our fits, we should investigate further before 
dismissing the SIS results in favour of the GIS results. 

Contamination of the GIS data 

The serendipitous source (hereafter assumed to be 
NSV 2321) seen in the SIS field is included in the CI S spe c- 
trum due to the larger PSF of the latter. In section 3.1.1 it 



21 Sept 1997 Observation 
. d.o.t. contour mop of o single-temp RS tit 



was remarked that NSV 2321 was intrinsically harder than 
the Iota Orionis spectrum (see also Fig. hi). To examine the 
effect of NSV 2321 on the fit results, a 1.06 arcmin radius 
circular region centred on NSV 2321 was extracted from the 
SISO periastron data. Although the quality of the spectrum 
was poor the fit results confirm that it is a harder source 
than Iota Orionis, which is consistent with the GIS fits hav- 
ing higher characteristic temperatures than the SIS fits. 

Intrinsic bias 

Another possibility for the differences between the SIS and 
GIS results may lie in the fact that compared with the 
SIS instruments, the GIS instruments are more sensitive at 
higher X-ray energies. This may result in fits to the GIS 




Figure 5. This figure shows the fit-statistic parameter space for a 
single-temp RS model fitted separately to the SISO and the GIS2 
periastron data. The contour levels plotted are the 90.0, 95.4, 
99.0, 99.73 and 99.99 per cent confidence regions. 



data intrinsically favouring higher characteristic tempera- 
tures. To investigate this we performed fits to a GIS spec- 
trum whilst ignoring the emission above selectively reduced 
energy thresholds. The results are shown in Fig. H where a 
clear bifurcation in the fit results can be seen. As more high 
energy X-rays are ignored, the fit jumps from a high temper- 
ature and low column to a low temperature and high column 
which is much more compatible with the SIS results. After 
discovering this behaviour the SISO data was analysed in a 
similar manner but reversing the process and selectively ig- 
noring X-rays below a certain threshold. Again, a bifurcation 
develops, this time towards the original GIS results. This is 
strong evidence that there is either an intrinsic bias in the 
instrument calibrations (either present in the original cal- 
ibrations, or a result of the instrument responses changing 
with time), or that the fitting process itself is fiawed i.e. that 
if we fit a single-temperature spectral model to an inherently 
multi-temperature source, the SIS fit will have a lower char- 
acteristic temperature than the GIS fit. Therefore are our 
fits simply telling us that the source is multi-temperature? 
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Table 4. Effective exposure times and background subtracted 
count rate and number of counts in the source region for the 
two ROSAT PSPC observations examined. The 'rf dataset was 
Boron filtered which blocks X-ray photons with energies between 
0.188-0.28 keV. 



Observation 



Phase 



Exp. 



Cts 



Ct rate 
(s-l) 



rp200700n00 
rf200700n00 



0.1697- 0.1841 
0.1841 - 0.1846 



746 
1305 



1153 1.545 
489 0.374 



In this case both the SIS and GIS results would be equally 
'incorrect'. 



3.2 Re-analysis of archival ROSAT data 

A search of the ROSAT public archive for observations of 
Iota Orionis yielded a total of 3 PSPC pointed observations. 
Two were back-to-back exposures with Iota Orionis on axis, 
of which one used the Boron filter. The other had Iota Orio- 
nis far off-axis (outside of the rigid circle of the window sup- 
port structure). An advantage of the ROSAT observations 
is that Iota Orionis is resolved from nearby X-ray sources. 
NSV 2321 is spatially resolved from Iota Orionis in two of 
the three PSPC pointings, and we have reanalysed both of 
these. Details of each observation are listed in Table M. 

A single temperature RS model could not reproduce the 
observed data for either observation. However, fits with two 
characteristic temperatures showed significant improvement 
and acceptable values of x? were obtained. Unfortunately an 
unrealistically soft component was obtained from the Boron 
filtered dataset which led to an artificially high luminosity. 
The fit to the other dataset is detailed in Fig. M and Table pi 
We note that if separate columns are fitted to each compo- 
nent the column to the harder component is much higher 
than that to the softer component and the estimated ISM 
value. 

Our results can be compared to the analysis by Geier 
et al. (1995), who obtained characteristic temperatures of 
0.12 and 0.82 keV from a two-temperature RS fit made to 
the rp200151n00 dataset. As already noted. Iota Orionis was 
substantially off-axis in this pointing and the nearby sources 
seen in Fig. Q were unresolved. This may account for the 
slightly higher temperatures that they obtained. 

We can also compare our 7? OSAT results to the ASCA 
results obtained in the previous section. In so doing we 
find distinct differences between the parameters fitted to 
the 3 individual detectors (SIS, GIS and PSPC), although 
roughly equivalent attenuated luminosities. In particular, 
the ROSAT spectrum has a significantly larger fiux at soft 
energies. 



3.3 The preferred datasets 

This leads us to the crucial question: which dataset do we 
have most confidence in? The negative points of each are: 

• The SIS datasets are invariably affected by the degra- 
dation of the CCDs. 
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Figure 6. This figure shows the fit-statistic parameter space for 
a single-temp RS model fitted to a GIS2 spectrum extracted from 
a 2 arcmin radius circle centered on Iota Orionis. The emission 
above certain energy thresholds was ignored during the fit process. 
For the top, middle and bottom panel these were 2.0, 1.7 and 
1.2 keV respectively. The contour levels plotted are the 68.3, 90.0 
and 99.0% confidence regions. The change from one minimum to 
the other is clearly shown. 
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Two temp RS fit 




ctiannel energy (keV) 

Figure 7. Two-temperature Raymond- Smith spectral model fit 
to the rp200700n00 ROSAT PSPC dataset. The large deviations 
between the model and the data at E ~ 0.5 keV may be due to 
uncertainties in the calibration of the PSPC detector response. 



• The GIS does not have a particularly good low energy 
response and does not spatially resolve the close surrounding 
sources which contaminate the subsequent analysis. 

• The PSPC has poor spectral resolution compared to 
the SIS and GIS. 

Despite the better spectral resolution of both instruments 
onboard ASCA, it would seem (from the uncertainties men- 
tioned in Section 3.1.4) that the ROSAT dataset returns 
the most accurate spectral fit parameters. However, the bet- 
ter spectral capabilities of ASCA convincingly demonstrates 
that there is no significant variation in the luminosity, spec- 
tral shape or absorption of the X-ray emission between the 
periastron and apastron observations. 



3.4 The X-ray luminosities 

The attenuated (i.e. absorbed) X-ray luminosity in the 0.5- 

2.5 keV band of the rp200700n00 ROSAT PSPC dataset 
(fitted by a two-temperature RS spectral model with one 
absorbing column) is ~ 1.6 x lO^^ergs"^. As previously 
mentioned it is often not easy to compare values with those 
estimated in previous papers because of the different spec- 
tral models assumed, the various energy bands over which 
the luminosities were integrated, the different methods of 
background subtraction employed, and the possible contam- 
ination from nearby sources (e.g. NSV 2321). However, as 
reported in Section bl Long & White (1980) estimated an at- 
tenuated luminosity of Lx = 2.3 x 10^^ ergs~^ from a single 
Einstein IPC observation. Chlebowski, Harnden & Sciortino 
(1989) obtained L^ = 3.0 x lO^^ergs"^ from seven Etn- 
stein IPC pointings. These values are within a factor of 2 
of the observed ROSAT 0.5-2.5 keV luminosity from the 
rp200700n00 dataset (1.6 x lO^^ergs"^). This is not un- 
expected given the differences in the various analyses. The 
0.5-10.0 keV luminosities from our analysis of the ASCA 
datasets (~ 1.0 x 10'^^ ergs~^ ), are also in rough agreement 
with the 0.5-2.5 keV 7? OSA rvalues, albeit slightly reduced 
for both the periastron and apastron observations. 

Because the ASCA datasets have much better spectral 



resolution, it would be a missed opportunity if we did not 
take advantage of this. Thus, whilst again acknowledging 
that the ASCA results may be systematically incorrect, in 
Fig. M we show the intrinsic and attenuated flux ratios of the 
periastron and apastron observations. The fluxes are inte- 
grated in the energy ranges 1.0-2.0, 2.0-3.0, 3.0-4.0, 4.0-5.0, 
5.0-6.0 and 6.0-7.0 keV although we note that the counts 
above 6.0 keV are minimal. Two-temperature RS spectral 
model fits were made to the combined SISO and SISl data, 
the combined GIS2 and GIS3 data, and all four datasets to- 
gether. The top panel shows the intrinsic flux ratio, whilst 
the lower panel shows the attenuated flux ratio. 

Immediately clear from Fig. W is the fact that the in- 
trinsic luminosity from the flt of all four datasets is almost 
constant between periastron and apastron, and this is also 
basically the case for the attenuated luminosity. If the ma- 
jority of the X-ray emission was due to colliding winds, one 
would expect the X-ray flux at periastron to be near maxi- 
mum, since Lx is proportional to 1/D. However, at energies 
above ~ 1 keV, the intrinsic Lx could be severely reduced 
around periastron as the shorter distance between the stars 
decreases the maximum pre-shock wind velocities. These ef- 
fects are not seen. 

Assuming that the bolometric luminosity of the 
Iota Orionis system is Li,oi — 2.5 x 10^ L0 (Stickland 
et al. 1987), we obtain log Lx/Lboi = —6.78 from the ob- 
served ROSAT luminosity. This compares to a value of 
— 6.601q j7 deduced from an Emstem observation (Chle- 
bowski, Harnden & Sciortino 1989). It is in even better 
agreement with the new results of Berghofer et al. (1997) 
who obtain log Lx/Ltoi ~ —6.81 (a — 0.38) for stars of lu- 
minosity classes III-V and colour index {B — V) < —0.25, 
which are appropriate for Iota Orionis. These results sug- 
gest that any emission from a colliding winds shock is at 
a low level (< 50 per cent of the total emission), consistent 
with the lack of variability between the ASCA pointings. We 
also repeat at this point the comment made by Waldron et 
al. (1998) on the validity of the Lx/Ltoi relationship: that 
the canonical ratio {Lx / LboisimlQ'"^) should only be inter- 
preted as an observed property of X-ray emission from OB 
stars. 



4 INTERPRETATION 

In the previous sections we have reported some very puzzling 
results. Iota Orionis was previously thought likely to have a 
colliding winds signature which would show significantly dif- 
ferent characteristics between periastron and apastron, and 
which would act as a test-bed for the latest colliding wind 
theories, such as radiative braking. Although there are con- 
cerns over the accuracy of the ASCA data and the absolute 
values of the fit results, it is still possible to compare rela- 
tive differences between the periastron and apastron obser- 
vations. However, we find that the ASCA X-ray lightcurves 
and fitted spectra show incredibly little difference between 
the two pointings. 

The extracted spectra also have lower characteristic 
temperatures than originally expected, especially for the 
apastron observation where the primary's wind should col- 
lide with either the secondary's wind or its photosphere at 
close to its terminal velocity (see Pittard 1998). 
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Figure 8. Ratio of the periastron and apastron intrinsic (top) and 
attenuated (bottom) fluxes as calculated from a two-temperature 
RS spectral model fltted to the ASCA instruments indicated. The 
symbols are plotted at the center of th e en ertfy band over which 
the fluxes were integrated (see Section p.4|). 



So, the issue is why there is not an obvious colliding 
winds X-ray signal in Iota Orionis. Given the evidence for 
colliding winds in other systems it seems unlikely that there 
is something fundamentally wrong with the colliding winds 
paradigm. Two interpretations are therefore possible. One is 
that the majority of the observed emission is from intrinsic 
shocks in the winds of one or both stars (c.f. Feldmeier et 
al. 1997a and references therein) which dominates a weaker 
colliding winds signal. In this scenario the variation of the 
observed emission with orbital separation would be expected 
to be small given the large emitting volume and the lack of 
substantial intrinsic variability. Alternatively the emission 
could be from a colliding winds origin but in a more com- 
plex form than one would naively expect. In the following 
subsections we focus our attention on each of these in turn. 



4.1 The circumstellar absorption 

Before discussing the merits and faults of each of these sce- 
narios, it is important to have a good understanding of the 
level and variability of the circumstellar absorption in the 
system. In Feldmeier et al. (1997a) a model for the wind at- 
tenuation of the primary star was derived from modelling 
of the UV/optical spectrum. They concluded that the wind 
of the primary 0-star is almost transparent to X-rays. If 
the attenuation is indeed negligible then (since the intrinsic 
X-ray colliding winds luminosity is expected to change as 
1/D) we would expect to see variability in the ASCA data. 



The fact that we do not has led us to re-consider the at- 
tenuation in the system from a colliding winds perspective, 
and we have calculated the absorbing column to the apex 
of the wind collision region as a function of orbital phase. 
This gives us an idea of the transparency of the wind along 
the line of sight to the shock apex. The basic assumptions 
involved in the construction of this model were: 

• The winds of both stars are spherically symmetric and 
are characterized by the velocity law 



1.0 - 







(1) 



where /3 = 0.8. 

• The wind collision shock apex is determined solely by a 
ram-pressure momentum balance. No radiation effects such 
as radiative inhibition or braking are included. 

• Where no ram-pressure balance exists the shock col- 
lapses onto the photosphere of the secondary and the apex 
is therefore located where the line of centers intersects the 
secondary surface. 

• The shock is not skewed by orbital motion. 

• The shock half-opening angle, 9, is determined from the 
equation 



V 



2/5 



.1/3 



for 10 ^ < J7 < 1 where 

M2V2 

V = —■ ■ 

M ivi 



(2) 



(3) 



Hi and V2 are the pre-shock velocities of the primary and 
secondary winds at the shock apex. This is basically the 
equation presented in Eichler & Usov (1993) but modified 
for non-terminal velocity winds. 

The line of sight to Earth from the shock apex is calculated 
and the circumstellar absorbing column along it evaluated 
from knowledge of the wind density and the position of the 
shock cone. The ISM absorption is not added because it is 
constant with orbital phase and we are primarily interested 
in the variation (and it is also much smaller at Nh ~ 2.0 x 
' cm-2). 

The orbit assumed for the calculation is shown in Fig. bl 
It is based on the cw_l model presented in Pittard (1998). 
The direction of Earth is marked, as well as the position of 
the secondary star relative to the primary at various phases. 
In Fig. hd we show the resulting orbital variability of the cir- 
cumstellar column with orbital and stellar parameters ap- 
propriate for Iota Orionis (the values used were again for 
model cw_l in Pittard 1998, which for convenience we list in 
Table pi) with a range of orbital inclinations. For i — 90° the 
numerical results were checked against an analytical form. 

Unlike binaries with circular orbits, the circumstellar 
column as a function of phase for an inclination i = 0° is 
not constant, being lower at apastron than at periastron. In 
Fig. [Kfl this variation is approximately an order of magni- 
tude. This is a direct result of the density of the primary's 
wind enveloping the shock cone being higher at periastron 
than at apastron, which mostly refiects the combination of 
the change in the orbital separation and the pre-shock ve- 
locity of the primary. 
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Figure 9. The assumed orbit in the frame of reference of the primary (marked with a cross at 0,0) for the cw_l model. The direction 
of Earth is marked with an arrow. At periastron the line of sight into the system may be through the secondary wind for a very short 
period. This is mainly dependent on the inclination angle of Iota Orionis, which is poorly known. The phases, (/>, over which the ASCA 
and i?05j4r exposures were made are delineated by diamonds. Also indicated are the phases corresponding to the images in Fig. 1 of 
Pittard (1998) and the lightcurve in Fig. ^. 



Table 5. Parameters used for the hydrodynamical models cw_l 
and cw_2 - see Pittard (1998) for a discussion of the observational 
estimates. A subscript 1 (2) indicates the value for the primary 
(secondary) star. 



Parameter 


Model cw_l 


Model cw_2 


Mi(M0yr-i) 


6.16 X lO-'' 


3.06 X lO-'^ 


M2(MQyr-i) 


1.88 X 10-« 


5.10 X 10-9 


!^cx=l(kms-i) 


2300 


2380 


'!'(x)2(kms~^ ) 


2200 


1990 



At higher inclinations (i > 30°) the shock apex is oc- 
culted by the secondary star just beforeperiastron. This is 
shown by a gap in the curves in Fig. hd. At i = 90° the 
primary star also occults the shock apex. On account of the 
large volume of any colliding winds X-ray emission these 
occultations are not expected to have any observational X- 
ray signatures. At inclinations i > 30° the maximum value 
of the circumstellar column becomes more asymmetric, and 
shifts in phase to where the shock apex is 'behind' the pri- 
mary (i.e. (/) ~ 0.15 - see Fig. bl). The orbital inclination has 
just been recalculated to lie in the range i = 50 — 70° (c.f. 
Marchenko et al. 2000). 

Comparing these results with the fitted columns in 
Table is difflcult, although they are of the correct or- 
der of magnitude. Assuming an ionized wind temperature 
of ~ 10* K, the T — 1 optical depth surface occurs at 
Nh ~ 2x 10^^ cm-2 and iVii ~ 6 x 10^^ cm'^ for 0.5 
and 1.0 keV photons respectively (c.f. for example, Krolick 
& Kallman 1984). This difference is mostly due to the Oxy- 
gen K-edge. Thus as Fig. p. 0| shows, the wind is indeed largely 
transparent at both of these energies, and for realistic values 
of the mass-loss rates. Only at phases near (j> = 0.0, where 



Nh ~ 10^^ cm-^, does the absorption begin to become ap- 
preciable (this is one of the causes of the X-ray minima 
in Figs. Illl and hsl). We therefore face two hard questions. 
Firstly, given the near transparency of the wind throughout 
most of the orbit, why don't we see the intrinsic variation 
in the colliding winds X-ray luminosity? Secondly, why do 
we not see enhanced absorption in our periastron observa- 
tion? In the hope of providing satisfactory answers to these 
questions we have calculated the expected colliding winds 
X-ray emission from a complicated hydrodynamical model, 
as detailed in the following section. 



4.2 The cw_l and cw_2 models 

In this section we investigate in greater depth whether a 
colliding winds model can be consistent with the lack of 
variability in the characteristic temperature and count rate 
of the ASCA data and the unremarkable X-ray luminosity 
of this system. In particular, we have two questions: i) is 
the X-ray luminosity from the wind collision low enough 
that the expected phase variability is lost in the intrinsic 
background; ii) are radiative braking effects stronger than 
anticipated and can they provide an explanation for both 
the low luminosity and the constant X-ray temperatures. 

Hydrodynamical models of the wind collision in 
Iota Orionis with ?7 = 31 and 72 (with terminal velocity 
values input as vi and V2) were presented in Pittard (1998). 
These simulations included the realistic driving of the winds 
by the radiation field of each star by using the line-force 
approximations of Castor, Abbott & Klein (1975). This ap- 
proach allows the dynamics of the winds and radiation fields 
to be explored, and a number of interesting and significant 
effects have been reported (e.g. Stevens & Pollock 1994; Gay- 
ley, Owocki & Cranmer 1997, 1999). We refer the reader to 
these papers for a fuller discussion. 

In Pittard (1998) it was found that the colliding winds 
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Orbital Phase 
Figure 10. Circumstellar column to the colliding winds shock 
apex as a function of orbital phase for the same stellar param- 
eters as used for model cw_l. The inclination of Iota Orionis is 
not well known so we show results for 4 sets of inclination (solid: 
i = 90°; dash: i = 60°; dot-dash: i = 30°; dot: i = 0°). The 
latest observational determination puts it at between 50 and 70° 
(Marchenko et al. 2000). Gaps in the curves indicate stellar oc- 
cultation. 



shock always collapsed down onto the surface of the sec- 
ondary star during periastron passage. On the other hand, 
there were large differences in the dynamics of the shock 
throughout the rest of the orbit. For model cw_2 the shock 
remained collapsed throughout the entire orbit, whilst for 
model cw_l the shock lifted off the surface of the secondary 
as the stars approached apastron, before recollapsing as the 
orbit headed back towards periastron. However, due to the 
necessary simplifications we cannot be sure that the shock 
actually behaves in this manner for a system correspond- 
ing to the particular parameters used (we also reiterate the 
interpretation of the latest optical data by Marchenko et 
al. (2000) who suggest that the shock appears to be detached 
even at = 0.0 — 0.1). Nevertheless, by comparing the or- 
bital variation of the X-ray emission from these models with 
the (lack of) observed variability by ASCA we may be able 
to determine if the observed emission is consistent with a 
colliding winds origin and if so whether the collision shock 
is actually detached from the secondary star near apastron. 
We refer the reader to Pittard & Stevens (1997) for spe- 
cific details of the X-ray calculation from the hydrodynamic 
models. We note that an orbital inclination of i = 47° was 
assumed for these calculations (c.f. Gies et al. 1996), and we 
do not expect them to change significantly in the light of 
the latest optical determination (i = 50 — 70°, Marchenko 
et al. 2000). 



4.3 An intrinsic wind shock interpretation 

Fig. [01 shows the attenuated luminosity from model cw_2 as 
a function of phase. Immediately apparent is its low over- 
all level in comparison to the .ROSAT measured luminosity 
{Lx ~ 2.0 X lO^'^ ergs~^ ). This gives us our first important 
insight into this system - the mass-loss rates of the two stars 
may be low enough with respect to other colliding wind sys- 
tems that the resulting emission is washed out against the 



intrinsic background. If this is indeed the case one would not 
even expect to notice the severe drop in the colliding winds 
emission predicted by the model around periastron. Simi- 
larly one would not expect to notice the variation in the 
characteristic temperature of the interaction region shown 
in Fig. Fl2. We conclude, therefore, that the colliding winds 
X-rays may be so weak that they are washed out by a rel- 
atively constant component from the intrinsic X-rays from 
each wind. However, we need to be sure that the emission 
measure of this intrinsic component (EMwind) is compara- 
ble to the X-ray observations (EMx-ray), since if we lower 
M , we lower EM^i„d and could run into the problem where 



EMu 



EM„ 



< EMx^ 



EM^ 



can be calculated from: 



47rr p {r)f{r)dr, 



(4) 



where /(r) accounts both for the volume filling factor and 
the differences in density of the shocked material relative to 
the ambient wind density p{r) (which assumes a smooth 
wind, given by p = M/47rr^ii). Previous wind instabil- 
ity simulations (Owocki et al. 1988; Owocki 1992; Feldmeier 
1995) are broadly consistent with a constant or slowly de- 
creasing filling factor as a function of radius (although de- 
tails of the wind dynamics are still largely unknown). Feld- 
meier et al. (1997a) further argue that due to the processes of 
shock merging and destruction a monotonically decreasing 
or roughly constant filling factor is appropriate. Hence, in 
the following a constant / will be supposed. For an instabil- 
ity generated shock model. X-ray emission is only expected 
once the wind has reached a substantial fraction (> 50 per 
cent) of Woo. We shall therefore take R,nin = 1.57?,. This is 
also sensible given that the observed X-ray flux doesn't show 
variability - if the emission occurred too close to the star 
eclipses would occur. With this value we also note that the 
fraction of the shocked wind occulted by the stellar disc is 
small. Feldmeier et al. (1997b) have also demonstrated that 
the X-ray emission from radii greater than 30i?, is negligi- 
ble, so we take this as our upper bound, Rmax- This is again 
reinforced by the fact that in Table H, Nh > Nism- Equn. H 
then becomes EM^i^d = f x EMsmooth, where EMsmooth is 
the emission measure of the smooth wind evaluated between 

^min anCl Umax • 

Derived filling factors for the shock-heated gas are of 
order 0.1-1.0 per cent for 0-stars (Hillier et al. 1993), but 
may possibly approach unity for near-MS B-stars (Cassinelli 
et al. 1994). For model cw_l, EMsmooth = 7.6 x 10'^'' cm"^ 
(assuming the velocity law in Equn. hi with (3 = 0.8). 
When compared to the values of EMx-ray in Table W we 
see that we require / ~ 0.001 — 0.42. For model cw_2, 
EMsmooth = a.Ox 10^* cm"^, and we require / ~ 0.02-1.07. 
Despite the large ranges these are compatible with the ob- 
servations (given that they also include the difference in den- 
sity between the shocked material and a smooth wind) , and 
therefore are consistent with an intrinsic wind shock inter- 
pretation. 

However, a possible problem with this interpretation 
arises from the larger flux seen in the ROSAT spectrum at 
soft energies with respect to the ASCA spectra. If we be- 
lieve that the spectral model parameters from the ROSAT 
and ASCA data are inconsistent, and that the variation in 
circumstellar attenuation is negligible, then this must be due 
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to a real change in the intrinsic flux. This appears somewhat 
untenable though, given the overall lack of variability in sin- 
gle 0-stars. 



4.4 A colliding winds interpretation 

In Fig. |l3| we show the X-ray lightcurve from model cw_l. 
This time the colliding wind emission is much stronger (ow- 
ing to the higher mass-loss rates assumed - see Table ph, 
and is again very variable. At the phases corresponding to 
the ASCA data the synthetic lightcurve is at local min- 
ima, and consistent with the observed luminosityrl The 
model lightcurve also predicts that the observed luminos- 
ity at phase — 0.17 should be higher than those at phases 
(f> — 0.0 and 0.5, with the latter at roughly the same levelPl 
From Table H we indeed find that this is the case, with 
the ROSAT luminosity matching the predicted lightcurve 
almost exactlyflJ. One slight worry is that the two Einstein 
observations at phases ~ 0.90 and 0.97 do not show a 
significant enhancement in count rate with respect to the 
observation at phase = 0.59. It would therefore be clearly 
useful to obtain new observations at these phases. 

The orbital variation of the characteristic temperature 
fitted to the colliding winds region of model cw_l is shown 
in Fig. [14. Our best estimate of the temperature varies from 
~ 0.2 keV at periastron to ~ 0.3 — 0.4 keV at apastron. This 
temperature variation is much less than expected from a 
simple consideration of the variation of the pre-shock veloc- 
ities along the line-of-centres, and reflects the fact that the 
fitted temperature is an average over the entire post-shock 
region. The temperatures predicted by the model are com- 
parable with the observed data and we thus conclude that 
this model demonstrates that the observed emission does 
not allow us to discard a colliding winds interpretation. 



4.5 Other possibilities/factors 

Another intriguing possibility is that radiative braking in the 
Iota Orionis system is actually quite efficient (moreso than 
the predictions of both models cw_l and cw_2). This might 
also explain the inference of Marchenko et al. (2000) that 
the shock is lifted off the surface of the secondary around 
periastron. This would require strong coupling between the 
B-star continuum and the primary wind, but cannot be ruled 
out by our present understanding. In the adiabatic limit 
(which is a good approximation for the wind collision in 
Iota Orionis at apastron), Stevens et al. (1992) determined 
that 



+ We note that the equations of Usov (1992) give un-attenuated 
luminosities roughly an order of magnitude below the results of 
our numerical ealeulations. 

8 There is some uncertainty in the true luminosity during the pe- 
riastron passage because the assumption of axisymmetry is poor 
at these phases, and additionally the optical data suggests that 
the wind collision region does not collapse onto the surface of the 
secondary. 

'I We note that the quoted GIS Lx {<j> = 0.5) is much higher 
than the GIS L^ {4> = 0.0) value. We conclude that this is due to 
uncertainties in the spectral fit parameters because the GIS back- 
ground subtracted count rates are nearly identical from Table H. 



Attenuated Luminosity (0.4-10.0 keV) 
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Figure 11. The 0.4—10.0 keV attenuated lightcurve for the col- 
liding wind model cw_2. A large drop in the emission around 
periastron is clearly seen. 



Lx (X M^v-''-^D-'^{l + ^f)/7f 



(5) 



where it was assumed that the temperature dependence of 
the cooling curve was A oc T~^'^ . This latter assumption 
is appropriate for post-shock gas in the temperature range 
~ 10^ — lO'^ K, which corresponds to pre-shock velocities of 
up to ~ 1000 km s~^ for solar abundance material. (Note: 
there is a typographical error in the corresponding equation 
of Stevens & Pollock, 1994). Assuming that r] is small (i.e. 
that the wind of the primary dominates) we then find that 



to first order, L^ 



Thus the X-ray luminosity is al- 



most inversely proportional to the pre-shock velocity of the 
primary wind. A reduction of 10 (50) per cent in the latter 
leads to an increase in the emission of 13 (130) per cent. 
Conversely, we find from the equations in Usov (1992), that 
Lx oc v~^''^ for the shocked primary wind and Lx oc u^ " for 
the shocked secondary wind. For model cw_l the emission 
from the latter appears to be dominant. Hence, for the same 
variations in the pre-shock velocity as discussed above, the 
emission is reduced by 5 and 30 per cent respectively. The 
discrepancy between these two sets of results is due to the 
different assumptions on the form of the cooling curve. Usov 
(1992) assumed cooling dominated by bremsstrahlung (i.e. 
A oc r^"), which is more appropriate at higher temper- 
atures (r > 10^ K). Clearly then, the situation is rather 
confusing at present, not least because the actual position 
of the colliding wind region relative to the two stars is very 
poorly known, as is the relative mass-loss rates of the stars. 
It is therefore difficult for us to be any more quantitative 
without performing a rigorous parameter-space study. 

It is also possible that the characteristics of the emis- 
sion may be altered by other physical processes. In investi- 
gating how electron thermal conduction may affect colliding 
wind X-ray emission, Myasnikov & Zhekov (1998) discov- 
ered that pre-heating zones in front of the shock have the 
overall effect of increasing the density of the wind inter- 
action region. The resulting X-ray emission was found to 
change markedly, with a large increase in luminosity and 
a significant softening of the spectrum. Softer X-rays suf- 
fer significantly higher absorption so the resulting observed 
emission may have an unremarkable luminosity. Whether 
the observed emission would show orbital variability is not 
clear at this stage. 

Finally, complex mutual interactions between the vari- 
ous physical processes occurring in colliding winds systems 
may also significantly alter the resultant emission. For in- 
stance, Folini & Walder (2000) suggest that the effects of 



© 2000 RAS, MNRAS 000, §-^ 



X-ray Analysis of the Eccentric 0-star Binary Iota Orionis 15 




0.4 0.6 

Phose 
Figure 12. The orbital variation of tiie cliaractcristic temper- 
ature from a one-temperature RS model fitted to the colliding 
wind model cw_2. Fits assuming solar abundance are marked with 
crosses, and those with a global floating abundance are marked 
with diamonds. A minimum in the characteristic temperature 
around periastron is clearly seen, reflecting the reduced pre-shock 
velocity of the primary wind in the hydrodynamical model. 
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Figure 13. The 0.4—10.0 keV attenuated lightcurve for the col- 
liding wind model cw_l. The small crosses (joined by the solid 
line) were calculated from the frames in Fig. 1 of Pittard (1998) 
at the phases indicated on Fig. & A large drop in the emission 
around periastron is clearly seen. However, when the shock de- 
taches from the secondary in this model a much larger grid is 
required to 'capture' all of the X-ray emission — this is shown 
by the dotted line. Also marked on this figure are the observed 
luminosities from the ASCA and i?05j4T pointings. The uncer- 
tainty in the count rates is a minimum of 1.2 per cent for the 
SISO periastron data, rising to a maximum of 4.3 per cent for 
the GIS2 apastron data. However, due to the uncertainties in the 
spectral shapes we conservatively assume a it 20 per cent error 
in their resultant luminosities. A very good fit is obtained if one 
retains a healthy scepticism about the exact size of the periastron 
minimum in the model, where our assumption of axisymmetry is 
least valid. We also mark on the relative count rates of the three 
Einstein observations reported by Snow, Cash &: Grady (1981). 
Because the conversion from count rate to luminosity is poorly 
known for these observations we have assumed uncertainties of ± 
35 per cent and arbitrarily adjusted their scaling to best fit our 
model curve. 



thermal conduction and radiative braking may positively 
reinforce each other. We are a long way from performing 
numerical simulations which would investigate this. 




Phose 
Figure 14. The orbital variation of the characteristic temper- 
ature from a one-temperature RS model fitted to the colliding 
wind model cw_l. Fits assuming solar abundance are marked with 
crosses, and those with a global floating abundance are marked 
with diamonds. A minimum in the characteristic temperature 
around periastron is clearly seen, reflecting the reduced pre-shock 
velocity of the primary wind in the hydrodynamical model. Note 
that unrealistically high values (> 1 keV) are obtained over the 
range 0.2 ^ <!> ^ 0.9 due to the grid not capturing all of the emis- 
sion. The rectangles and ellipses show the correct results of solar 
and floating abundance fits to a larger hydrodynamical grid over 
this phase range. Thus the characteristic temperature is within 
the range 0.2 - 0.4 keV over the entire orbital phase. 



5 COMPARISON WITH OTHER COLLIDING 
WIND SYSTEMS 

Most early-type systems with strong colliding wind sig- 
natures are WR-I-OB binaries. These are different from 
OB-I-OB systems in a number of ways which may explain 
why it appears that stronger colliding wind signatures are 
obtained from WR-I-OB systems. First, the high values of 
M wR provide more wind material which can be shocked. 
Second, the spectral emissivity for WC wind abundances is 
greater than that for solar abundances at the same mass- 
density (i.e. mass-loss rate - see Stevens et al. 1992). How- 
ever, both of these points may also act in reverse (high mass- 
loss rates also provide more absorption, and the emissivity 
of WN wind abundances for a given mass density is below 
that of solar). At this point it is therefore instructive to 
compare the X-ray characteristics of Iota Orionis with other 
early-type binary systems. The latter can be subdivided into 
three distinct groups as detailed in the following subsections. 



5.1 Those with strong coHiding winds signatures 

Into this group fall the well-known binaries WR 140 
(HD 193793) and 7^ Velorum (WR 11, HD 68273). Both 
these systems show clear evidence for colliding stellar winds 
including phase-variable emission and a hard spectrum. 
WR 140 (WC7 -I- 04-5, P=7.94 yr) is famous for its episodic 
dust formation during periastron passage (Williams 1990). 
Strong X-ray emission from WR 140 was discovered by EX- 
OSAT and its progressive extinction with phase by the WC7 
wind has been used to derive the CNO abundances of the 
WR wind (Williams et al. 1990). Possible non-thermal X-ray 
emission has also recently been discovered in this system, as 
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witnessed by the variability of the Fe K line (Pollock, Cor- 
coran & Stevens 1999). 

ROSAT observations of 7^ Velorum (WC8 + 09I, 
P=78.5 d) were presented by Willis, Schild & Stevens 
(1995). The hard X-ray flux showed a phase repeatable in- 
crease by a factor of 4 when the system was viewed through 
the less-dense wind of the O-star, which was attributed to 
the addition of a harder component to the largely unvarying 
softer emission. More recent ASCA observations (Stevens et 
al. 1996, and a re-analysis by Rauw et al. 2000) confirmed 
these results. A detailed comparison of this data with syn- 
thetic spectra generated from a grid of hydrodynamical col- 
liding wind models provided evidence that the hard X-ray 
emission comes directly from the wind collision. 

5.2 Those showing some evidence of colliding 
winds 

The best example of this category is V444 Cyg (WR 139, 
HD 193576), a well-studied eclipsing WN5 -I- 06 binary with 
an orbital period P—4.21 d and well known physical param- 
eters. Corcoran et al. (1996) confirmed from phase-resolved 
-ROSylr observations that although the luminosity was 1-2 
orders of magnitude lower than the predictions of Stevens et 
al. (1992), there was an orbital dependence of the flux. This 
on its own did not completely rule out the O-star being the 
source of the X-rays, though fairly convincing evidence for 
colliding winds remained. First, the high value of the char- 
acteristic temperature outside of eclipse {kT ~ 1.7 keV) 
was much greater than the typical temperatures of single 
O-stars {kT ~ 0.5 keV). Second, the variation of the emis- 
sion was consistent with a shock location near the O-star 
surface, as would be expected given the higher mass-loss 
rate of the WR-star and previous deductions from ultra- 
violet /[/£ observations (Shore & Brown 1988). Finally, the 
assumption of instantaneous terminal wind velocity overes- 
timates the observed X-ray luminosity. More recent ASCA 
(Maeda et al. 1999) and optical (Marchenko et al. 1997) ob- 
servations of this system have reinforced this interpretation. 
However, some puzzling X-ray characteristics remain such 
as large amplitude short time-scale variability which is not 
easily explained by a colliding winds model (Corcoran et 
al. 1996). 

Another system which shows some, but not overwhelm- 
ing, characteristics of colliding winds is 29 UW Canis Ma- 
joris (HR 2781; HD 57060). This is also a short-period 
(P=4.3934 d) eclipsing binary system. ROSAT ^-my obser- 
vations presented by Berghofer & Schmitt (1995b) showed 
phase-locked variability with a single broad trough centered 
near secondary minimum. However a colliding winds in- 
terpretation was dismissed by the authors on the assump- 
tion that the wind momenta of the component stars were 
broadly similar and thus that the shocked region was nearly 
planar (Wiggs & Gies 1993). Such a scenario should form 
a double-peaked lightcurve (e.g. Pittard & Stevens 1997) 
and the majority of the emission was therefore attributed 
to the secondary's wind. The unremarkable X-ray luminos- 
ity also seemed in accordance with single O-star emission. 
However, if the wind momenta are much more imbalanced 
than believed, the lightcurve would be expected to have a 
single broad minimum, in agreement with the observations. 
HD 57060 also suffers from the Struve-Sahade effect (Stick- 



land 1997; Gies, Bagnuolo & Penny 1997), and whilst the 
cause of this remains uncertain, colliding stellar winds is one 
possible interpretation. We therefore conclude that colliding 
winds may still be signiflcant in this system. 

Other strong candidates for colliding winds emission in- 
clude HD 93205 and HD 152248 (c.f. Corcoran 1996), both 
of which are due to be observed with the XMM satellite, and 
HD 165052 (see Pittard & Stevens 1997 for more details). 

5.3 Those with no colliding winds signature 

Many early-type binaries observed by X-ray satellites show 
no signature of colliding winds emission. For the most part 
this is because the observations are generally short, the 
lightcurves extremely sparse and the expected luminosi- 
ties near or below previous detection limits. However, one 
fairly bright frequently observed system which shows no 
signs whatsoever of colliding winds emission is S Orionis A 
(HD 36486), a spectroscopic and eclipsing binary (09.5 II -f- 
BO HI) with a 5.7 day period. As far as we know there has 
been no published paper which specifically investigates any 
potential colliding winds emission in this system. 



6 CONCLUSIONS 

In this paper we have analysed two ASCA X-ray observa- 
tions of the highly eccentric early-type binary Iota Orio- 
nis, which were taken half an orbit apart at periastron and 
apastron. Although all previous observations were short and 
without good phase coverage, it was expected that a strong 
colliding winds signature would be seen. This expectation 
was based on the almost order of magnitude variation in 
the orbital separation between these phases and the known 
dependence of the X-ray luminosity of colliding winds (and 
suspected dependence of the spectral shape) with this. In 
turn it was hoped that this would reveal further insight 
into the physics of radiative driving and stellar winds from 
early-type stars, and in particular the additional interac- 
tions possible in binary systems. Indeed, Iota Orionis was 
selected as a target precisely because substantial changes 
were expected which would allow us to 'probe' the dynam- 
ics of the wind collision and infer the amount of radiative 
braking/ inhibition. 

Although our analysis was complicated by various prob- 
lems experienced with the ASCA datasets, we found the 
emission to be surprisingly constant and unvarying between 
the two observations, in direct contradiction with our expec- 
tations. A further analysis of archival ROSAT data, demon- 
strated a relatively low X-ray luminosity from this system 
of Lx ~ 1.0 — 1.6 X 10^^ ergs~^ . Using a simple model we 
confirmed that the column along the line of sight to the ex- 
pected position of the shock apex was unlikely to produce 
any significant attenuation of the shock X-rays, except per- 
haps around periastron. Based on an expected 1/D varia- 
tion in the emission this did not help to explain our constant 
count rate, and subsequently significantly more complex hy- 
drodynamical models were used to investigate the problem. 

Despite strongly phase- variable emission from the mod- 
els, both were consistent with the observations. The model 
with the lower mass-loss rates (cw_2) predicted attenuated 
luminosities an order of magnitude below the observational 
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Table 6. Summary of the strengths and weaknesses of the two interpretations of the X-ray emission considered. 



Intrinsic Wind Shocks 



CoUiding Winds 



Pros Lack of variabiUty in ASCA data. 



Lightcurve and spectra of model cw_l are consistent 
with the observational data, despite neither matching 
our initial simple expectations. 



Low mass-loss rates of the stars can reduce CW emission 
but without running into problems with the X-ray emission 
measure exceeding that of the wind. 



Cons Circumstellar attenuation around periastron may be 
significant. 



Modelling is very difficult and contains a number of 
assumptions and approximations. 



Optical data suggests that the wind collision shock 
does not crash onto the surface of the secondary 
during periastron passage. 



data, implying that intrinsic wind shocks emit the majority 
of the observed emission. On the other hand, the model with 
the higher mass-loss rates (cw_l) predicted minima in the 
colliding winds emission which coincided with the phases of 
the ASCA observations, and slightly enhanced emission at 
the phase of the ROSAT observation. This was also con- 
sistent with the data, implying that the observed emission 
could also be interpreted as purely colliding winds emission. 
In Table m we summarize the strengths and weaknesses of 
each of these. 



Unfortunately it is impossible to distinguish between 
these two interpretations with the limited dataset currently 
available. Additional observations are clearly needed. New 
data with increased spatial resolution (to avoid source con- 
fusion from objects such as NSV 2321), spectral resolution 
and counts (to constrain spectral and luminosity variabil- 
ity), and phase coverage is necessary if we are to begin dis- 
tinguishing the competing interpretations. Observations at 
other wavelengths to further refine the system parameters 
(particularly the mass-loss rates of both stars) are also de- 
sirable. 



Finally we compare the Iota Orionis observations with 
those obtained for other early-type binary systems. Al- 
though Iota Orionis is not the first well-studied early-type 
binary to show no clear signature of colliding winds emission 
{5 Orionis A is another), it is the first binary with a highly 
eccentric orbit not to do so. The lack of a clear colliding 
winds signature in Iota Orionis makes it noticeably differ- 
ent from other eccentric binaries (e.g. WR 140, 7^^ Velorum) 
where such signatures are readily seen at X-ray energies. 
Thus, Iota Orionis raises new open questions in our under- 
standing of early-type binaries. In particular, we would like 
to know why the colliding winds emission is not obvious in 
Iota Orionis. Is it simply a result of the low mass-loss rates, 
a particular combination of system parameters, or are there 
deeper forces at work? This may be the first indication of 
the importance of radiative braking on a line-driven wind. 
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